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ABSTRACT 

Context. In the last five years the Fermi Large Area Telescope (LAT) instrument detected GeV y-ray emission from five novae. The 
GeV emission can be interpreted in terms of an inverse Compton process of electrons accelerated in a shock. In this case it is expected 
that protons in the same conditions can be accelerated to much higher energies. Consequently they may produce a second component 
in the y-ray spectrum at TeV energies. 

Aims. We aim to explore the very-high-energy domain to search for y-ray emission above 50 GeV and to shed light on the acceleration 
process of leptons and hadrons in nova explosions. 

Methods. We have performed observations with the MAGIC telescopes of the classical nova V339 Del shortly after the 2013 outburst, 
triggered by optical and subsequent GeV y-ray detections. We also briefly report on VHE observations of the symbiotic nova YY Her 
and the dwarf nova ASASSN-13ax. We complement the TeV MAGIC observations with the analysis of contemporaneous Fermi -LAT 
data of the sources. The TeV and GeV observations are compared in order to evaluate the acceleration parameters for leptons and 
hadrons. 

Results. No significant TeV emission was found from the studied sources. We computed upper limits on the spectrum and night-by- 
night flux. The combined GeV and TeV observations of V339 Del limit the ratio of proton to electron luminosities to L p < 0.15 L e . 

Key words. Novae: cataclysmic variables - Gamma-rays: stars - Binaries: general - Stars: activity 


1. Introduction 

A classical nova is a thermonuclear runaway leading to the ex¬ 
plosive ejection of the envelope accreted onto a white dwarf 
(WD) in a binary system in which the companion is either fill¬ 
ing or nearly filling its Roche surface (fBode & Eva~nsl 120081 : 
IStarrfield et alj20l2tlWoudt & Ribeirol2014l) . They are a type of 


* Corresponding authors: J. Sitarek (jsitarek@uni.lodz.pl), 
W. Bednarek (bednar@uni.lodz.pl), R. Lopez-Coto (rlopez@ifae.es) 


cataclysmic variable, i.e. optically variable binary systems with 
a mass transfer from a companion star to WD. Novae are typi¬ 
cally detected first in optical observations when the brightness of 
the object increases by 7-16 magnitudes. A thermal X-ray con¬ 
tinuum is often seen in the energy spectra of novae. Symbiotic 
novae, like classical novae, are initiated by a thermonuclear ex¬ 
plosion on the surface of the WD. However in the case of symbi¬ 
otic novae the WD is deeply imme rsed in the wind of a late-type 
companion star (see e.g. lShore et al.ll20TTll2012 ). 
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The diffusive shock acceleration at the blast wave of sym¬ 
biotic novae was expected to accelerate part icles up to ener¬ 
gies of a few TeV dTatischeff & Hernanzll2007T) . In 2010 the first 
GeV y-ray emission was detected by the Fermi Large Area Tele¬ 
scope (LAT) from the symbiotic nova V407 Cyg (lAbdo et al.l 
l2010t) . The y-ray emission can be explained in terms of lep- 
tonic or hadronic models (lAbdo et al.l 1201 Ot lAckermann et al.l 
120141) . The local radiation fields create a target for the inverse 
Compton (IC) scattering of the electrons. Protons accelerated 
in the same conditions can interact with the matter produc¬ 
ing y-rays via proton-proton interactions. Several models have 
been put forward to expla in the GeV radiation. For instance, 
ISitarek & Bednarekl (12012 1 attribute the GeV y-ray emission to 
the IC process on the strong radiation field of the red giant. The 
same model predicts a second component in the TeV range due 
to proton-p roton interactions with the wind of the red giant. Al- 
ternativelv. lMartin & Dubusl (120131) consider acceleration of lep¬ 
tons and hadrons in the nova shock. In that model the magnetic 
field, which determines the acceleration efficiency, is obtained 
assuming an equipartition with the thermal energy density up¬ 
stream of the shock resulting in the maximum energy of protons 
estimated to be ~ 300 GeV. The GeV y-ray emission is then ex¬ 
pected mostly from leptonic processes, namely the IC scattering 
of the nova light by the electrons accelerated in the shock. 

In the last three years Fermi-LAT has discovered GeV 
y-ray emission from four additional no vae, VI32 4 Sco, 
V959 Mon, V339 Del and VI369 Cen dCheunget al.ll2013bl: 
lAckermann et al.l l2014ll . In addition Fermi-LAT has reported 
V745 S co and Nova Sgr 2015 No. 2 as lower significance can¬ 
didates tlCheung et al.ll2014al 120151) . Most of these sources are 
classical novae. Contrary to the symbiotic ones, they do not ex¬ 
hibit a strong wind of the companion star, but still they all show 
similar spectral properties. In classical novae particle accelera¬ 
tion can occur, for example, in a bow shock between the nova 
ejecta and the interstellar medium or in weaker internal shocks 
due t o inhomogeneity of the nova ejecta (lAckermann et alJ 
120141 ). In particular, the orbital motion of the system can shape 
the nova ejecta into a faster polar wind and a slower region 
of de n ser material al o ng the equatorial plane (IChomiuk et alJ 
12014 . Metzger et ah (2015 ) suggest that the y-ray emission 
might come from hadronic interactions of a faster outflow with 
an ejected shell. 

So far no very-high-energy (VHE; £ >100 GeV) y-ray emis¬ 
sion has been detected from any nova event. VERITAS observa¬ 
tions of the symbiotic nova V407 Cyg beginning 10 days after 
the nova explosion yielded a differential upper li mit o n the flux 
at 1.6 TeV of 2.3 x 10 12 erg cnr 2 s 1 2 (lAliu et al.ll20H) . 

Since late 2012 the MAGIC collaboration has been conduct¬ 
ing a nova follow-up program in order to detect a possible VHE 
y-ray component. At first the program focused on symbiotic no¬ 
vae. After the reports of the detection of GeV emission from 
classical novae by the Fermi -LAT the program was extended 
also to bright classical novae and opened up to additional out¬ 
bursts from other cataclysmic variables. 

In this paper we report on the observations performed with 
the MAGIC telescopes of V339 Del and present an updated 
analysis of Fermi-LAT data contemporaneous with those obser¬ 
vations. In Section [2] we describe the MAGIC telescopes and 
the Fermi -LAT instrument. The observations of V339 Del are 
presented in Section [3 In Section [4] we discuss the results of 
the GeV-TeV observations of V339 Del in terms of a hadronic - 
leptonic model. In addition to the classical nova V339 Del, 
MAGIC observed also dwarf nova ASASSN-13ax and symbi¬ 
otic nova YY Her. Neither GeV nor TeV y-ray emission was 


detected from those two sources. The upper limits on these two 
objects are summarized in the online-only Appendix[A]while the 
main text of the paper focuses on V339 Del. 


2. Instruments 

V339 Del and the other outbursts observed by MAGIC were first 
detected and observed by optical instruments. The analysis of 
quasi-simultaneous Fermi -LAT observations provides additional 
context for the MAGIC results. 

2.1. MAGIC telescopes 

The VHE y-ray observations were obtained using the MAGIC 
telescopes. MAGIC is a system of two 17 m Cherenkov tele¬ 
scopes located on the Canary Island of La Palma at a height 
of 2200 m a.s.l. (lAleksic et al.li2016al) . The telescopes record 
y-rays with energies above ~50GeV. The sensitivity of the 
MAGIC telescopes in the best energy range (> 300 Ge V) is ~ 
0. 6% o f Crab Nebula flux in 50 h of observations (lAleksic et al.l 
l2016bl) . The data were analyzed using the standard analysis 
chain (IZanin et al .1120131: lAleksic et al.ll2016bl ). The significance 
of a y-ra y excess was computed according to Eq. 17 from 
iLi & Mai (119831) . The upp er limits on the flu x were calculated 
following the approach of iRolke et al.l (120051) using 95% confi¬ 
dence level (C.L.) and accounting for a possible 3 0% systematic 
uncertainty on the effective area of the instrument (lAleksic et al.l 

I20l6bh . 

2.2. Fermi-LAT 

The Fermi- LAT, a space-based, pair-conversion telescope, de- 
tects photons with energies from 20MeV to > 300 GeV 

dAtwood et al.l l2009h . We analyzed the LAT data in the energy 
range lOOMeV - 300 GeV using an unbinned maximum likeli¬ 
hood method dMattox et al.ll 19961) as implemented in the Fermi 
Science Tools v9r32p5. We applied the P7REP_SOURCE_V15 
LAT Instrument Response Functions (IRFs) and used the as¬ 
sociated standard Galactic and isotropic diffuse emission mod¬ 
els matched to the Pass 7 reprocessed Source class event se- 
lectioifl We selected events within a region of interest (ROI) 
of 15° in radius cen t ered o n the LAT best-fit position reported 
bv lAckermann et ah] (2014) for V339 Del and required a maxi¬ 
mum zenith angle of 100° to avoid contamination by Earth limb 
photons. Because some of the LAT data were acquired during 
pointed mode observations, we applied an appropriate filteiQ, se¬ 
lecting good quality data at times when either the rocking angle 
was less than 52° or the edge of the analysis region did not ex¬ 
ceed the m aximum z enith angle at 100°. Sources from the 2FGL 
catalogue (iNolan et al.ll2012» located within the ROI were in¬ 
cluded in the model used to perform the fitting procedure. 


3. Observations of V339 Del and results 

V339 Del was a fast, classical CO nova detected by optical 
observations on 2013 August 16 (CBET #3628), MJD 56520. 
The nova was exceptionally bright reaching a magnitude of 


1 The P7REP data, IRFs, and diffuse models (gll_iem_v05_revl .fit and 
iso_source_v05.txt) are available at http://fermi.gsfc.nasa.gov/ssc 

2 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
Cicerone_Likelihood/Exposure.html 
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Fig. 1 . Multi-wavelength light curve of V339 Del during the outburst 
in August 2013. Top panel: Optical observations in the V band obtained 
from the AAVSO-LCC0 service. Middle panel: The Fermi- LAT flux 
(filled symbols) and upper limits (empty symbols) above 100 MeV in 1- 
day (circles, thin lines) or 3-day (squares, thick lines) bins. A 95% C.L. 
flux upper limit is shown for time bins with TS<4. Bottom panel: Upper 
limit on the flux above 300 GeV observed with MAGIC telescopes. The 
gray band shows the observation nights with MAGIC. The dashed gray 
line shows a MAGIC observation night affected by bad weather. 

a http://www.aavso.org/lcg 


V~ 5 mag (see top panel of Fig. Q}, and it triggered follow¬ 
up observations at frequencies ranging from radio to VHE y- 
rays. Photometric measu reme n ts su ggest a distance for V339 
Del of 4.5 + 0.6 kpc (ISchaefer et al.ll2014l) . Nearly a month after 
the optical detection. X-ray emission was de tecte d in the 1-10 
keV energy band by Swift/XRT ( Page et al . 2013 ). Afterwards, 
the object became a low energy X-ray source with most of the 
photons detected in the 0.3-1 keV energy range (lOsborne et alJ 
120131 ). The object shows large amplitude variations and a 54 s 
quasi-periodic oscillation in the soft X-ray energy band. These 
are possibly explained by the spin of the white dwarf or an 
oscillation in th e nuclear burning rate dBeardmore et al I [20 1 3b 
iNess et al]l2013h . The spectroscopic observations performed on 
2013 August 18 revealed emission wings extending to about 
±2000 km/s and a Balme r absorptio n component at a velocity 
of (600 + 50) km/s dShore et al.ll2Q13bb . The pre-outburst optical 
images revealed the progenitor of nova V339 Del to be a blue 
star ( Denisenko et alJl2013l) . 


Originally, MAGIC observations of V339 Del were moti¬ 
vated by its extreme optical outburst. The subsequent detec¬ 
tion of GeV emission by the Fermi -LAT from the direction of 
V339 Del added incentive for VHE observations. MAGIC ac¬ 
quired data starting on the night of 2013 August 16 but these 
were marred by poor weather. The good quality data used for 
most of the analysis spans 8 nights between 2013 August 25 
and September 4. The total effective time was 11.6 h. In addition 
to the nightly upper limits we performed a dedicated analysis 
of the poor-quality (affected by Calima, a dust layer originating 
from Sahara) night of 2013 August 16. We applied an estimated 
energy and afilection area correction based on LIDAR measure¬ 
ments dFruck et al.1120141) . No VHE y -ray signal was found from 
the direction of V339 Del. We computed night-by-night integral 
upper limits above 300 GeV (see bottom panel of Fig. Q]) and dif¬ 
ferential upper limits for the whole good quality data set in bins 
of energy (see Section[4]). 

Nova V339 Del was the subject o f a Fer mi Target of Op¬ 
portunity (ToO) observation (lHavs et al.ll2013l) triggered by the 
optical discovery (CBET #3628); a pointed observation favoring 
the nova started on 2013 August 16 and lasted for 6 days. The 
y-ray emission from V339 Del wa s first detected by Fermi -LAH 
in a 1-day bin on 2013 August 18 dAckermann et al.ll2014l) . The 
emission peaked on 2013 August 22 and entered a slow decay 
phase afterwards (Fig. |T]). We fitted the flux for the light curves 
shown in the middle panel of Fig. |T|by assuming a power-law 
spectral model with the normalization left free to vary and the 
photon index fixed to a single value. We selected a fixed value 
for the photon index of 2.3 by calculating the average of the most 
significant 1-day bins (Test Statistic values TS>9]j. 

The spectral energy distribution (SED) for V339 Del shown 
in Fig. [2] was extracted in 5 logarithmically spaced energy bins 
from 100 MeV to 100 GeV. Similarly to the light curves, energy- 
binned data shown in Fig. [2] were fitted using a power law and 
calculating a 95% C.L. upper limit for bins with TS<9. In the 
period coincident with the MAGIC observations (2013 August 
25 to September 4), the Fenni-LXT spectrum can be described 
by a power law with an index of 2.4 + 0.2 and flux above 100 
MeV of (0.15 ± 0.04) x 10~ 6 phcm~ 2 s -1 . Only statistical error 
is included since that dominated over systematic error in these 
data. The spectral fit for this period had a TS of 49 and did 
not permit a constraint on an exponential cut-off in energy. The 
Fermi-LAT analysis for the full decay phase, 2013 August 22 
to September 12 (MJD 56526-56547), provided a more signifi¬ 
cant signal with a TS of 121 and a similar value of flux above 
100 MeV, (0.13 ± 0.03) x 10 _6 phcm _2 s _1 . The spectrum for 
the longer time period was fitted by an exponentially cut-off 
power law with an index of 1.4 ± 0.3 and a cut-off energy of 
1.6 + 0.8 GeV. The fit improved in significance by 3.3cr with re¬ 
spect to a power-law model. The most energetic photon associ¬ 
ated with V339 Del (89% probability for the best-fit model) had 
an energy E — 5.9 GeV and was recorded on 2013 August 30, 
i.e., within the time period covered by the MAGIC observations. 

4. Modeling of y- rays from nova V339 Del 

Of the 3 objects observed by the MAGIC telescopes and dis¬ 
cussed in this paper, V339 Del is the only one detected at GeV 
energies by the Fermi- LAT. Extensive optical observations pro¬ 
vided constraints on the companion star and the photosphere of 
the nova. Therefore it has the highest potential for constraining 

3 The source significance (in sigmas) is ~ yjTS assuming one degree 
of freedom 
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Fig. 2. Differential upper limits on the flux from V339 Del as measured by MAGIC (filled squares) and the flux measured by Fermi -LAT (empty 
crosses) in the same time period, 2013 August 25 to September 4. The thin solid line shows the IC scattering of thermal photons in the nova’s 
photosphere. The dashed line shows the y-rays coming from the decay of n° from hadronic interactions of the relativistic protons with the nova 
ejecta. The dotted line shows the contribution of y-rays coming from IC of e + e~ originating from n + n~ decays. Thick solid lines show the total 
predicted spectrum. The total energy of electrons is 6 x 10 41 erg and the assumed proton to electron luminosity ratio is L p IL e = 0.1. Electrons and 
protons are inj ected w ith a power law with a spectral index of 1.5 and the cut-offs reported in the figures. Photosphere parameters (see Table QJ 
obtained from lMetzger et al ; 12015 ?). 


the leptonic and hadronic processes in novae, and we concen- 
trate the modeling efforts o n it. We follow the modified model of 
ISitarek & Bednarek](12012!) . The original scenario was applied to 
the symbiotic nova V407 Cyg. In that case the GeV y-ray emis¬ 
sion was attributed to IC of the electrons on the strong radiation 
field in the vicinity of the red giant companion star. For V339 
Del, however, the radiation field of the companion star is not as 
strong. The photosphere of the nova provides a dominant target 
for the IC process. Moreover, the wind of the companion star is 
not as dense as in the V407 Cyg symbiotic system. Nevertheless, 
if protons of sufficient energy are accelerated in the nova shock, 
they can interact with the ejecta of the nova produ cing pions that 
decay to y-rays. As in lSitarek & Bednarek ( 2012 ) the GeV y-ray 
emission can be used to constrain the parameters describing the 
acceleration of the electrons, which are otherwise poorly known. 
Protons will then be accelerated in the same conditions but up to 
much higher energies due to lower energy losses. 

In order to apply the above model we need to evaluate the 
radiation field which is encountered by electrons at the time of 
the observations by MAGIC and Fermi- LAT. First, we need to 
estimate the parameters of the nova photosphere about 10 days 
after the optical detection. The spectral evolution of the V339 
Del was studied bv lMetzger et all ( 2015 ). During the time of the 
observations by the MAGIC telescope the reported temperature 
of the photosphere was ~ 7000 K. In the early phases of the nova 
(fireball stage) the optical observations were consistent with a 
pseudophotosphere with a temperature of about 10000K. With 
a reported optical luminosity of ~ 6 x 1O 4 L 0 this corresponds 
to a radius of the photosphere of 1 .2 x 10 13 cm (see Table |T}. 
Following the approach of 1 Shore et alJ (l2013al) we also used a 
second set of photosphere parameters, in which UV emission, 


not available at this late stage of V339 outburst is estimated using 
measurements of a similar source, namely OS And. 

Let us consider a photosphere of a radius R p h — 
10 I3 i? P A,i3 cm and a temperature T p i, = 10 4 7' 4 K. The leptons 
and hadrons are accelerated at a distance of R s h - 10 li R s h,n cm 
from the photosphere. Assuming that the velocity of the shock 
is ~ 1000 km/ s, similar to the one observed in another GeV- 
emitting nova (IChomiuk et al.l20l4l) . we estimate that during the 
MAGIC observations the distance from the WD to the shock was 
R s h ~ 10 I4 cm. The y-rays with energy E y = 10.EioGeV can be 
produced via IC scattering in the Thomson regime of thermal 
photons by electrons with energy of 

E e - 22x E 1 { q 2 / (77t(l + cos/I)) 1/2 [GeV], (1) 

where fi is the angle between the electron and the direction to the 
point on the photosphere where the thermal photon was emitted. 
On the other hand, by comparing the energy losses from the IC 
scattering with the accelerati on rate we obtain a maxim um en¬ 
ergy for electrons of (see e.g. ISitarek & BednarekH2012l) 

Ee, max = l3(^ 4 B) ll2 R sh /(T 2 4 R ph ) [GeV], (2) 

where B is the magnetic field at the shock (measured in Gauss), 

and £ ; = 10 4 f. 4 is the acceleration coefficient. The acceleration 

coefficient is defined by the acceleration time r acc = 1 E/£ . 4 fi. 

By comparing the above two formulae we obtain 

UB = 2.9E w TlR 1 ph l{R) h {\ + cos/I)). (3) 

In the same conditions protons with energy E p (measured in 
units of GeV) can be accelerated with a time scale of 

EpR 2 , (1 + cos B ) 

Tacc, p = 1 V«Ufl) = 0-34 £ . , [s], (4) 

£ 10 T 4 R ph 
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Table 1 . Parameters characterizing the optical emission of V339 Del 
(photosphere temperature T\ x 10 4 K, radius R p h,ij, x 10 13 cm and lumi¬ 
nosity L) according to the two scenarios assumed in the modeling of the 
GeV and TeV emission. 



r 4 

R P h,n 

L/L b 

Metzger et al. (2015) 

0.7 

1.2 

6x 10 4 

optical+UV 

1.3 

0.4 

8 x 10 4 


The acceleration can be limited by the dynamic time scale of 
tj ~ 10 days (after which time most of the MAGIC observations 
were performed) or by the energy losses of protons from the pp 
collisions. The time scale of the latter can be computed as 

t pp = ( cr pp n H kc ) -1 , (5) 

where cr pp « 3 x 10~ 26 cm 2 is the interaction cross section, is 
the density of the nova ejecta, the inelasticity coefficient, k ~ 0.5, 
is the fraction of energy lost in each interaction, and c is the 
speed of light. The density of the ejecta will decrease as the nova 
shock progresses with the speed of v = 10 3 t ; 3 km s _1 following 

n H = 4.4 x 10 12 AU 5 /(v 3 f 3 ) [cm- 3 ], (6) 

where 10 5 M 5 M (: , is the total mass ejected during the outburst. 
Thus, 

Tpp = 500(v 3 f^)/M_5 = 780tf 3 /u3 /M^5 [s], (7) 


where we used that R s h = v(86400 G) « 8.6 x 10 I3 cm. 

By comparing the acceleration and cooling time scales we 
obtain 


E p = 2300 


EioT 4 R 2 phl3 R s h,i3 

(1 + COSj8)M_5 


[GeV]. 


( 8 ) 


Usi ng the para meters of the photosphere obtained by 
(Metzger et alJ (2075!) (see TableQJ a fortnight after the peak and 
taking into account the break in the GeV spectrum at ~ 1.6 GeV, 
it is plausible to expect protons to be accelerated at least up to 
energies of ~ 1.6 TeV. The maximum energy of the protons is ex¬ 
pected to be lower, 1.1 TeV, in the case of the second set of pho¬ 
tosphere parameters computed from the combined optical and 
UV fit. 

The acceleration can be also limited by the dynamic time 
scale instead of the hadronic interaction losses. By comparing 
the dynamic time scale with the proton cooling time scale we 
find that the former starts to dominate only after 13A7 7 2 v] 3/2 
days. Therefore the accelerated protons will mostly cool down 
on the time scale of the MAGIC observations due to energy 
losses in hadronic interactions. The normalization of both com¬ 
ponents is determined by L p /L e , i.e. the ratio of the total power 
of accelerated protons to that of elec trons. _ 

We used the numerical code of ISitarek & Bednareki (12012l) 
to model the Fermi -LAT GeV spectrum and to compare the sub- 
TeV predictions with the MAGIC observations. We consider that 
the electrons and protons accelerated in Fermi-like acceleration 
obtain a power-law energy spectrum with a spectral index of 1.5. 
The spectra of electrons and protons cut off at energies deter¬ 
mined by equations [2] and [8] respectively. In Fig.[2]we show the 
predictions for leptonic or hadronic spectra compared with the 
Fermi-L AT and MAGIC measurements. 

The Fermi-L AT spectrum can be described mostly by IC 
scattering of the thermal photons in the nova’s photosphere by 


electrons. The expected hadronic component overpredicts the 
MAGIC observations at ~ 100 GeV by a factor of a few for the 
case of equal power of accelerated protons and electrons (i.e. 
L p = L,,). Using the upper limits from the MAGIC observa¬ 
tions we can place the limit on L p < 0.15L e . We checked that 
our results are unchanged when using radiation field parameters 
resulting in a significantly different emission in the UV range, 
which is not constrained by the observations of the V339 nova at 
the time of the MAGIC observations. 

The increased power in electrons compared to protons may 
be related to how particles with different mass are injected in the 
acceleration process. Interestingly, the appearance of energetic 
e + e- pairs from nuclear decays produced in the nova explosion 
could help to inject them preferentially into the shock acceler¬ 
ation process. On the other hand. [Schlickeiserl ( 2003 ) suggests 
that in a low-beta plasma acceleration of electrons is preferred 
over protons if the particles are accelerated out of a thermal pop¬ 
ulation. Both effects could lower the L p !L e ratio. 


5. Conclusions 

No VHE y-ray emission was found from the direction of V339 
Del. The contemporaneous Fermi-LAT observations revealed 
GeV emission from V339 Del. We modeled this GeV emission 
as the IC of thermal photons from the photosphere by the GeV 
electrons accelerated in the nova shock. We used the Fermi -LAT 
and MAGIC observations of V339 Del to constrain the number 
of protons accelerated in the same conditions as the electrons 
in the nova shock. The modeling shows that the total power of 
accelerated protons must be < 15% of the total power of acceler¬ 
ated electrons. MAGIC will continue to observe promising y-ray 
nova candidates in the following years. 

Acknowledgements. The MAGIC Collaboration would like to thank the Insti¬ 
tute de Astroffsica de Canarias for the excellent working conditions at the Ob- 
servatorio del Roque de los Muchachos in La Palma. The financial support 
of the German BMBF and MPG, the Italian INFN and INAF, the Swiss Na¬ 
tional Fund SNF, the ERDF under the Spanish MINECO, and the Japanese 
JSPS and MEXT is gratefully acknowledged. This work was also supported 
by the Centro de Excelencia Severo Ochoa SEV-2012-0234, CPAN CSD2007- 
00042, and MultiDark CSD2009-00064 projects of the Spanish Consolider- 
Ingenio 2010 programme, by grant 268740 of the Academy of Finland, by 
the Croatian Science Foundation (HrZZ) Project 09/176 and the University 
of Rijeka Project 13.12.1.3.02, by the DFG Collaborative Research Centers 
SFB823/C4 and SFB876/C3, and by the Polish MNiSzW grant 745/N-HESS- 
MAGIC/2010/0. and NCN 2011/01/B/ST9/00411. JS is supported by Fundacja 
UL. We thank M.A. Perez-Torres for the information about nova ASASSN-13ax. 
The Fermi LAT Collaboration acknowledges generous ongoing support from a 
number of agencies and institutes that have supported both the development and 
the operation of the LAT as well as scientific data analysis. These include the 
National Aeronautics and Space Administration and the Department of Energy 
in the United States, the Commissariat a l’Energie Atomique and the Centre Na¬ 
tional de la Recherche Scientifique / Institut National de Physique Nucleaire et de 
Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto 
Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT), High Energy Accelerator Research Organiza¬ 
tion (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the 
K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish 
National Space Board in Sweden. Additional support for science analysis during 
the operations phase is gratefully acknowledged from the Istituto Nazionale di 
Astrofisica in Italy and the Centre National d’Etudes Spatiales in France. Au¬ 
thors would like to thank S. N. Shore for scientific discussions and providing an 
optical spectrum of V339 Del and UV spectrum of OS And. We would also like 
to thank the anonymous referee for comments that helped to improve the paper. 


References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, Science, 329, 817 
Ackermann, M., Ajello, M., Albert, A., et al. 2014, Science, 345, 554-558 


Article number, page 5 of[7] 

















A&A proofs: manuscript no. magicnova 


Aleksic, J., Ansoldi, S., Antonelli, L. A., et al. 2016a, Astroparticle Physics, 72, 
61-75 

Aleksic, J., Ansoldi, S., Antonelli, L. A., et al. 2016b, Astroparticle Physics, 72, 
76-94 

Aliu, E., Archambault, S., Arlen, T., et al. 2012, ApJ, 754, 77 
Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 2009, ApJ, 697, 1071 

Bode, M. F., & Evans, A. 2008, Classical Novae, 2nd Edition. Edited by 
M.F. Bode and A. Evans. Cambridge Astrophysics Series, No. 43, Cam¬ 
bridge: Cambridge University Press, 2008. 

Beardmore, A. P., Osborne, J. P., Page, K. L 2013, The Astronomer’s Telegram 
#5573 

Castro-Tirado, A. J., Martin-Carrillo, A., & Hanlon, L. 2013, The Astronomer’s 
Telegram, #5314 

Cheung, C. C., Jean, P, & Shore, S. N. 2013, The Astronomer’s Telegram, # 5653 
Cheung, C. C., Jean, P., & Shore, S. N. 2014, The Astronomer’s Telegram, # 5879 

Cheung, C. C., Jean, P., & Shore, S. N., 2015, The Astronomer’s Telegram, 
#7283 

Chomiuk, L., Linford, J. D., Yang, J., et al. 2014, Nature, 514, 339 

Copperwheat, C. M., Bersier, D. F., Shappee, B. J., et al. 2013, The Astronomer’s 
Telegram, #5195 

Deacon, N. R„ Hoard, D. W., Magnier, E. A., et al. 2014, A&A, 563, A129 
Denisenko, D., Jacques, C., Pimentel, E., et al. 2013, IAU Circ., 9258, 2 

Fruck, C., Gaug, M., Zanin, R., et al. 2014, Proc of 33rd ICRC, Rio de Janeiro, 
Brazil, arXiv: 1403.3591 

Hays, E., Cheung, T., Ciprini, S., The Fermi LAT Collaboration 2013, The As¬ 
tronomer’s Telegram # 5302 

Herbig, G. H. 1950, PASP, 62, 211 

Li, T.-R, & Ma, Y.-Q. 1983, ApJ, 272, 317 

Martin, P„ & Dubus, G. 2013, A&A, 551, A37 

Metzger, B. D„ Finzell, T„ Vurm, I„ et al. 2015, MNRAS, 450, 2739 

Mattox, J. R., Bertsch, D. L., Chiang, J., et al. 1996, ApJ, 461, 396 

Munari, U., Rejkuba, M.. Hazen, M., et al. 1997, A&A, 323, 113 

Munari, U., Dallaporta, S., Castellani, F., et al. 2013, The Astronomer’s Telegram 
#4996 

Ness, J. U., Schwarz, G. J., Page, K. L., et al. 2013, The Astronomer’s Telegram 
#5626 

Nolan, P. L., Abdo, A. A., Ackermann, M.. et al. 2012, Astrophys. J. Supp. 199, 
31 

Osborne, J. P., Page, K., Beardmore, A. et al. 2013, The Astronomer’s Telegram 
#5505 

Page, K. L. and Beardmore, A. P. 2013, The Astronomer’s Telegram # 5429 
Rolke, W. A„ Lopez, A. M„ & Conrad, J. 2005, NIM A, 551, 493 
Schaefer, G. H., Brummelaar, T. T., Gies, D. R., et al. 2014, Nature, 515, 234 

Schlickeiser, R. 2003, Energy Conversion and Particle Acceleration in the Solar 
Corona, 612, 230 

Shore, S. N., Wahlgren, G. M., Augusteijn. T., et al. 2011, A&A, 527, AA98 
Shore, S. N., Wahlgren, G. M., Augusteijn. T., et al. 2012, A&A, 540, AA55 

Shore, S. N., Schwarz, G. J., Alton, K., et al. 2013a, The Astronomer’s Telegram, 
#5409 

Shore, S. N., Skoda, R, Korcakova, D., et al. 2013b, The Astronomer’s Telegram, 
#5312 

Shore, S. N. 2013, The Astronomer’s Telegram, # 5410 
Sitarek, J., & Bednarek, W. 2012, Phys. Rev. D, 86, 063011 

Stanek, K. Z., Shappee, B. J., Kochanek, C. S., et al. 2013, The Astronomer’s 
Telegram, #5186 

Stanfield, S., Iliadis, C., Timmes, F. X., et al. 2012, Bulletin of the Astronomical 
Society of India, 40, 419 

Tatischeff, V., & Hemanz, M. 2007, ApJ, 663, L101 

Woudt, P. A. & Ribiero, V. A. R. M. (eds) 2014, Stella Novae: Past and Future 
Decades. ASP Conference Series, Vol. 490 (Astron. Soc. Pacific) 

Zanin, R., Carmona, E., Sitarek, J., et al., 2013, Proc of 33rd ICRC, Rio de 
Janeiro, Brazil, Id. 773 


1 ETH Zurich, CH-8093 Zurich, Switzerland 

2 Universita di Udine, and INFN Trieste, 1-33100 Udine, Italy 

3 INAF National Institute for Astrophysics, 1-00136 Rome, Italy 

4 Universita di Siena, and INFN Pisa, 1-53100 Siena, Italy 

5 Croatian MAGIC Consortium, Rudjer Boskovic Institute, University 
of Rijeka and University of Split, HR-10000 Zagreb, Croatia 

6 Saha Institute of Nuclear Physics, 1\AF Bidhannagar, Salt Lake, 
Sector-1, Kolkata 700064, India 

7 Max-Planck-Institut fur Physik, D-80805 Miinchen, Germany 

8 Universidad Complutense, E-28040 Madrid, Spain 

9 Inst, de Astroflsica de Canarias, E-38200 La Laguna, Tenerife, 
Spain 

10 University of Lodz, PL-90236 Lodz, Poland 

11 Deutsches Elektronen-Synchrotron (DESY), D-15738 Zeuthen, 
Germany 

12 IFAE, Campus UAB. E-08193 Bellaterra, Spain 

13 Universitat Wurzburg, D-97074 Wurzburg, Germany 

14 Centro de Investigaciones Energeticas, Medioambientales y Tec- 
nologicas, E-28040 Madrid, Spain 

15 Universita di Padova and INFN, 1-35131 Padova, Italy 

16 Institute of Space Sciences, E-08193 Barcelona, Spain 

17 Technische Universitat Dortmund, D-44221 Dortmund, Germany 

18 Unitat de Fisica de les Radiacions, Departament de Fisica, and 
CERES-IEEC, Universitat Autonoma de Barcelona, E-08193 Bel¬ 
laterra, Spain 

19 Universitat de Barcelona, ICC, IEEC-UB, E-08028 Barcelona, 
Spain 

20 Japanese MAGIC Consortium, ICRR, The University of Tokyo, De¬ 
partment of Physics and Hakubi Center, Kyoto University, Tokai 
University, The University of Tokushima, KEK, Japan 

21 Finnish MAGIC Consortium, Tuorla Observatory, University of 
Turku and Department of Physics, University of Oulu, Finland 

22 Inst, for Nucl. Research and Nucl. Energy, BG-1784 Sofia, Bulgaria 

23 Universita di Pisa, and INFN Pisa, 1-56126 Pisa, Italy 

24 ICREA and Institute of Space Sciences, E-08193 Barcelona, Spain 

25 Universita dell’Insubria and INFN Milano Bicocca, Como, 1-22100 
Como, Italy 

26 now at Centro Brasileiro de Pesquisas FIsicas (CBPF\MCTI), R. Dr. 
Xavier Sigaud, 150 - Urea, Rio de Janeiro - RJ, 22290-180, Brazil 

27 now at NASA Goddard Space Flight Center, Greenbelt, MD 20771, 
USA and Department of Physics and Department of Astronomy, 
University of Maryland, College Park, MD 20742, USA 

28 Humboldt University of Berlin, Institut fur Physik Newtonstr. 15, 
12489 Berlin Germany 

29 now at Ecole polytechnique federale de Lausanne (EPFL), Lau¬ 
sanne, Switzerland 

30 now at Finnish Centre for Astronomy with ESO (FINCA), Turku, 
Finland 

31 also at INAF-Trieste 

32 also at ISDC - Science Data Center for Astrophysics, 1290, Versoix 
(Geneva) 

33 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 


Article number, page 6 of[7] 



M. L. Ahnen et al.: Very-high-energy y-ray observations of novae and dwarf novae with the MAGIC telescopes 


Appendix A: (Online only) YY Her and 
ASASSN-13ax 

Two additional objects, the symbiotic nova YY Her and dwarf 
nova ASASSN-13ax, were observed as part of the broader cam¬ 
paign but neither was detected to have GeV or TeV emission. In 
this appendix we report the values of the differential upper limits 
from those two sources. The Fermi-L AT analysis was config¬ 
ured as described for V339 Del with the ROIs centered on the 
optical positions for YY Her and ASASSN-13ax. Data were se¬ 
lected within the time windows reported, preceding the optical 
peak and containing approximately the following two weeks of 
observations. The region model was derived from the 2FGL cat¬ 
alogue and contains all sources within 20° of the ROI center. The 
spectra for YY Her and ASASSN-13ax were described using a 
power law with the photon index fixed to 2.2. The flux normal¬ 
ization for the source of interest and the normalization for the 
isotropic template were left free to vary while the other source 
parameters were fixed, either to catalogue values or refitted val¬ 
ues if the source flux was detected in the analysis window and 
deviated by more than lcr from the catalogue value. The upper 
limits were calculated at 95% C.L. using the Bayesian method 
provided with the Fermi Science ToolsQ. 

The upper limits from MAGIC observations were computed 
in 5 bins per decade in energy. The photon spectral index as¬ 
sumed for both sources was 2.6. 

YY Her is a symbiotic nova system that undergoes a re¬ 
current pattern of outbursts. MAGIC observations of YY Her 
occurred on the night of April 22, 7 days after the optical 
maximum. The observations were motivated by the symbiotic 
nature of the source. In the absence of detectable emission, 
flux upper limits at 95% C.L. were calculated to be 5.0 x 
10-“ phcm~ 2 s _1 above 300 GeV. Emission was not detected 
in the LAT over the interval 2013 April 10 to April 30 (MJD 

56392.5 to 56412.5). Upper limits at 95% C.L. were calculated 
to be 2.8 x 10~ 8 phcm~ 2 s _1 above 100 MeV. Differential upper 
limits obtained from the Fermi -LAT and MAGIC observations 
of Y Y Her are summarized in Tab. IA.1I 

ASASSN-13ax is a member of a different class of cata¬ 
clysmic variable, the dwarf novae, which are known for signifi¬ 
cantly weaker optical outbursts (typically 2-8 magnitudes) than 
classical novae. The source was observed by MA GIC d ue to 
its very strong optical outburst of 7.7 magnitudes (IStanek et al.l 
l2013l f. Instead of undergoing a thermonuclear explosion on the 
surface of the WD, these outbursts are caused by the gravita¬ 
tional energy release from an instability in the accretion disk 
surrounding the WD. The MAGIC observations were performed 
on two consecutive nights starting on 2013 July 4, soon after the 
optical outburst seen on 2013 July 1. In the absence of detectable 
VHE emission, upper limits at 95% C.L. were calculated to be 

1.5 x 10~ 12 phcuT 2 s -1 above 300 GeV. Emission was not de¬ 
tected in the LAT over the interval 2013 June 25 to July 15 (MJD 

56468.5 to 56488.5). Upper limits at 95% C.L. were calculated 
to be 1.6 x 10 -8 phcm~ 2 s~* above 100 MeV. Differential upper 
limits obtained from the Fermi-LA I and MAGIC observations 
of ASASSN-13ax are summarized in Tab. IA.2l 


Table A.l. Differential upper limits on the flux from YY Her as mea¬ 
sured by the Fermz'-LAT and MAGIC. The bins extend from E mln to E max 
while the upper limit value is computed at the energy of E UL , 



Fermi- LAT, MJD 

: 56392.5-56412.5 

Emin [GeV] 

Emax [GeV] Et 

jl [GeV] 

F ul [TeV cnrV 1 ] 

0.100 

0.316 

0.178 

6 .1x10-“ 

0.316 

1.00 

0.562 

4.2 xl0~ 12 

1.00 

3.16 

1.78 

4.7 xl0~“ 

3.16 

10.0 

5.62 

1 .2x10-“ 

10.0 

100.0 

31.6 

1.5x10-“ 

MAGIC, MJD=56405 

50.0 

79.2 

59.4 

7.1x10-“ 

79.2 

125.6 

94.2 

1.9x10-“ 

125.6 

199.1 

149.3 

3.8 xlO- 12 

199.1 

315.4 

236.6 

6.5xl0- 12 

315.4 

500.0 

375.1 

4.3xl0- 12 


Table A.2. Differential upper limits on the flux from ASASSN-13ax as 
measured by the Fermz-LAT and MAGIC (see text for details). Columns 
as in Table lA.il 



Fermi -LAT , MJD 56468.5 - 

56488.5 

Emin [GeV] 

Emax [GeV] 

Eul [GeV] 

Ful [TeV cm^s -1 ] 

0.100 

0.316 

0.178 

3.3x10-“ 

0.316 

1.00 

0.562 

3.7 xlO- 12 

1.00 

3.16 

1.78 

2.6 xlO- 12 

3.16 

10.0 

5.62 

1.4x10-“ 

10.0 

100.0 

31.6 

1.5x10-“ 

MAGIC, MJD: 56478-56479 

79.2 

125.6 

94.2 

1 .2x10-“ 

125.6 

199.1 

149.3 

1.2 xlO- 12 

199.1 

315.4 

236.6 

1.7xl0- 12 

315.4 

500.0 

375.1 

7.3x10-“ 

500.0 

792.4 

594.4 

l.lxlO- 12 

792.4 

1255.9 

942.1 

7.9x10-“ 

1255.9 

1990.5 

1493.1 

9.1x10-“ 

1990.5 

3154.7 

2366.4 

1.7x10-“ 

3154.7 

5000.0 

3750.5 

9.3x10-“ 


4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/python_tutorial.html 
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